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Abstract 

It has been recently argued that localized, unstable, but extremely long-lived 
configurations, called oscillons, could affect the dynamics of a first order elec¬ 
troweak phase transition in an appreciable way. Treating the amplitude and the 
size of subcritical bubbles as statistical degrees of freedom, we show that thermal 
fluctuations are not strong enough to generate subcritical configurations able to 
settle into a an oscillon long-lived regime. 
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Nontopological soliton solutions of classical field theories were introduced a number 
of years ago by Rosen Q] and their properties have been studied by many authors @]. 

Unlike magnetic monopoles and cosmic strings, which arise in theories with non¬ 
trivial vacuum topology, nontopological solitons are rendered stable by the existence 
of a conserved Noether charge carried by helds conhned to a hnite region of space. 
The minimum charge of the stable soliton depends upon ratios of coupling constants 
and in principle can be very small. 

Scenarios for actually producing such objects in the early Universe have also been 
discussed [^. In particular, this issue has been recently rekindled in ref. [Q where lo¬ 
calized, time-dependent, spherically-symmetric solutions of nonlinear scalar field the¬ 
ories, called oscillons, were studied and shown to be, although unstable, extremely 
long-lived. Indeed, their lifetimes can be of order of (10^ — 10^) m~^, where m is the 
mass of the scalar held, i.e. much longer than that for a conhguration obeying the 
Klein-Gordon equation for a free scalar held, of order of 5 m~^. 

Oscillons naturally appear during the collapse of spherically symmetric held con- 
hgurations: if a bubble is formed at rest at the time t = 0 with, for example, a 
’’Gaussian” shape 

0(x,O) = ae-l"l'/^', (1) 

i. e. with an amplitude a at its core and initial radius R, after having radiated most of 
its initial energy, the bubble settles into a quite long-lived regime, before disappearing 
by quickly radiating away its remaining energy. 

The conditions required for the existence of the oscillons are, apart from having the 
initial energy above a plateau energy, essentially two |Q: i) the initial amplitude of the 
held at the core needs to be above the inhection point of the potential in order to probe 
the nonlinearities of the theory and ii) the conhguration must have an initial radius R 
bounded from below. To explain these conditions fairly analitically (conclusions are 
conhrmed numerically), we can follow ref. [Q and consider the potential 
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A solution 0(x, t) to the equation of motion has energy 
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Since during the oscillon regime the subcritical conhguration is characterized by a 
slowly varying radius, we can model the oscillon by writing 
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where the radius R^ is kept hxed (a good approximation supported by numerical 
analysis). 





The equation of motion for a{t) can be linearized writing a{t) = a{t) + 6a{t), so 
that the fluctuation 6a{t) satishes the linearized equation 

Sa = —uj'^{a,R^:)6a, 

2/- T-) \ 3a/2 2 4\/6 - 2/1 3 \ , , 

uj^{a,R^) = Xa^ - — aoma + \1 +—^ j . (5) 

For a;^(a, i?*) < 0 fluctuations about a are unstable, driving the amplitude away from 
its vacuum value and thus avoiding the rapid shrinking of the initial conhguration. 
These are the fluctuations responsible for the appearance of the oscillon and its relative 
long lifetime [Q. 

For hxed ao, uj^(a,R^) does have a minimum for hmin — 0.51 {aom/X), hence 
oscillons are possible only for 

R^ > i?min — y^ 3/ (0.28ao - A) m. (6) 

For R^ > u;^(a, R^) will be negative for amplitudes 

a- < a < a+, (7) 


where 


Sv^aom , \/2 

a± = -;— ± — 

27 A 3 


96 O-Qm? 




1/2 


( 8 ) 


81 A2 

In the limit of very large i?*, i?* S> 1/m, a_ becomes independent from R^. To give a 
numerical example, in the degenerate case ao = (3/a/2) A^^^ 
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where we have indicated with Oinf and Umax the inflection point closest to a = 0 and 
the maximum of the potential V{a), respectively. 

It is then clear why oscillons can form only if their initial amplitude at the core is 
above the inflection point Omf and why their initial radius R cannot be too small in 
order to feel the nonlinearities of the potential. 

One of the motivations for studying the evolution of unstable spherically-symmetric 
conhgurations comes from the original papers analyzing the role subcritical bubbles 
may play in the dynamics of weak hrst order phase transitions . Considering models 
with double-well potentials in which the system starts localized on one minimum, for 
sufficiently weak transitions correlation-volumes bubbles of the other phase could be 
thermally nucleated, giving rise to an effective phase mixing between the two available 
phases before the reaching from above of the tunneling temperature at which critical 
bubbles are expected to be nucleated. This could have dramatic consequences for any 
electroweak baryogenesis mechanisms |p. 











Although the presence of thermal fluctuations in any hot system is undisputed, 
their role in the dynamics of weakly first order phase transitions is still under debate 0]. 
However, it is clear that, i/thermally nucleated, long-lived oscillons could appreciably 
affect the dynamics of a weak first order phase transition at the weak scale. Although 
their lifetime is small in comparison with the expansion time-scale for temperatures 
T ~ 100 GeV, if oscillons are produced in large enough numbers, their presence will 
substantially increase the equilibrium number-density of subcritical bubbles of the 
broken phase. This could effectively make the transition weaker than what predicted 
from the effective potential. Also, instabilities on the expanding critical bubble walls 
could be generated by collisions with oscillons, implying that the usual assumption of 
spherical evolution of the walls may be incorrect. 

The aim of this Letter is to investigate whether oscillons can be really present at 
the onset of a first order electroweak phase transition, i.e if subcritical bubbles with 
initial amplitude at their core and initial radius R satisfying the above conditions i) 
and a) can be thermally nucleated and affect the usual picture of the phase transition 
dynamics. 

To answer this question we will treat both the initial amplitude and the size of 
subcritical bubble as statistical degrees of freedom along the same lines of what done 
by Enqvist et al. in ref. [0. 

First order phase transition and bubble dynamics in the Standard Model have 
lately been studied in much detail, and it has become increasingly clear that for 
Higgs masses considerably heavier than 60 GeV, the electroweak phase transition is 
only of weakly first order. For Higgs mass Mh > 100 GeV the calculations, both 
perturbative and lattice ones, confront technical problems and it is conceivable that 
for such large Higgs masses the electroweak phase transition is close to a second order 
and does not proceed by critical bubble formation. Therefore, in this Letter we use 
a phenomenological Higgs potential for the order parameter 0 suitable for a simple 
description of a first order phase transition: 

V(0, T) = \m\T)ct? - (10) 

The properties of the oscillons for the potential ( 0 ) can be easily derived from the 
analysis made for the potential (|^) with the substitution a^m aT. Namely, the 
oscillon stage can be obtained only if subcritical configurations have initial amplitude 
greater than the inflection point 0inf and sufficiently large size. We also expect that, 
when increasing A, the minimum necessary value for R increases, whereas the smallest 
available value of the amplitude at the core decreases [||. 

When discussing oscillons one has to be sure that initial configurations, which 
eventually will give rise to an oscillon, are not critical bubbles. Indeed, for the potential 
( pTl) and at tunneling temperature Tf, critical bubbles become solutions of the equation 


of motion: if they are nucleated with an initial radius Rc (or larger) they can grow 
converting the metastable phase 0 = 0 into the stable phase with lower energy. 

Most of the dynamical properties of the electroweak phase transition associated 
with the potential Eq. (0), such as the smallness of the latent heat, the bubble nucle- 
ation rate and the size of critical bubbles, have been discussed extensively in 0. For 
the purposes of the present paper it suffices to recall only some of the results. 

Assuming that there is only little supercooling, as seems to be the case for the 
electroweak phase transition, the bounce action can be written as 

a 

= A3/2 35 (A -1)2 ’ 


where A(T) = 9Xm‘^{T)/{2a‘^T^). The cosmological transition temperature is deter¬ 
mined from the relation that the Hubble rate equals the transition rate oc 
yielding S/Tf ~ ln(Mp[/Ty) ~ 150, where Tf is the transition temperature. Thus we 
obtain from Eq. ( 0 ) 

A(U) ~ 1 - 0.0442^ = 1 - A. (12) 

On the other hand, small supercooling implies that 1 — A = 5 -C 1, i.e. a -C 500A^/^. 
Solving for A in Eq. (|^) yields the transition temperature Tf. One hnds 

2(R 

m\Tf) = ^\{Tf)T] . (13) 

The extrema of the potential are given by 

0min,max(T) = (l ± • (14) 


Expanding the potential at the broken minimum 0min(E) we hnd 

1 1 / \ 

- t = Tf = - -Xft,,, = -0.00218 + O , (15) 

The height of the barrier is situated at 0max — 4>min/2. with V (0min,Ec) = 14iax = 
a'^T^^/(144 A^), where Tc is the temperature at which H(0) = H(0min), given by the 
condition = (2 T^/9 A). As ~ Tf we may conclude that the thin wall 

approximation is valid if —e/14iax = 0.314 <C 1, or a -C lOA^/^. Thus the 

small supercooling limit is clearly satished if the thin wall approximation is valid. 
The closest inflection point to 0 = 0 at T ~ is given by 

0inf ~ 0.42 0max- (16) 


To get the size of the critical bubble we still need the surface tension. One easily 
hnds 


a = d0^2 V{T,) 
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We define the critical bubble radius by extremizing the bounce action. The result is 
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Therefore Rc is much larger than the correlation length ^{Tf) = l/m{Tf) at the 
transition temperature, as it should. 

Let us hrst make the general observation that it is the actual transition temperature 
Tf rather than the critical temperature which is relevant for the study of oscillons. 
This is true in the sense that if oscillons are not important at Tf, they most certainly 
will not be so at Tc- As we shall show, it actually turns out that oscillons are not 
present even at Tj. This justihes, in retrospect, our choice T = Tf for performing the 
calculations. 

In the case of a weak hrst order phase transition the critical bubble is typically well 
described by a thin wall approximation, where the conhguration has a hat ’highland’ 
(with 0 determined by the non-zero minimum of the potential) and a steep slope 
down to 0 = 0. Therefore it seems natural that also a large subcritical bubble should 
resemble the critical one, i.e. when R increases, the form of the subcritical bubble 
should deform smoothly so that, when R = R^, the bubble is a critical one. 

Motivated by this observation, let us dehne a subcritical bubble as a functional of 
both the amplitude a and the radius R. For this purpose one has hrst to study the 
behaviour of the potential as a function of the amplitude. At Tj there is a interval 
0 G [a_, a+] where IL(0) < 0. If the amplitude of the bubble is in that interval, 
there exists a critical bubble-solution of the bounce action. This means that we have 
a relation Rc = Rc{a) which reproduces Eq. (|T8|) if a = 0min- Therefore Rc{a) serves 
as an upper limit for the initial radius i? of a subcritical bubble in that region: if 
R > Rc{a) we exclude such a conhguration since it should give rise to a critical bubble 
and not, eventually, to an oscillon with hnite lifetime. 

These considerations lead us to dehne diherent Ansatze for various regions in the 
(a, i?)-plane. When 0 G [a_, a+], we use an Ansatz such that when R —> Rc{a), 
the held conhguration goes towards the thin wall form. For small R we use a simple 
gaussian conhguration. For other values of a we always take a thin-wall like Ansatz. 
Thus we write for 0 G [a_, a+] and R < Rc{a) 



(19) 


where t is the time coordinate^ and 



^Note, however, that we need not to specify the explicit time evolution of a and R when dealing 
with statistical averages. 
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Such an Ansatz reproduces the requirement that when R —> Rc{a), subcritical bubbles 
should resemble critical ones. In practise the statistical averages depend only weakly 
on a because the main contribution to them comes from the region of small a amd 
large R. Therefore we assume for simplicity that criticality depends weakly on a and 
take Rc{a) = Rc to be a constant whenever possible. 

For (f) ^ [a_, a+] we assume that no gaussian component is present and write 
simply 

R) = a(t)0t(i?). (23) 


However, the statistical averages are expected to be quite insensitive of the precise 
form of the conhguration. 

These Ansatze can be plugged into the action 


S[a, R] = J 


X 




VW 


(24) 


from which the Lagrangian in terms of the dynamical variables a and R can be ex¬ 
tracted. In the practical calculation we have, whenever possible, approximated by 
the step function. After that is a simple matter to calculate the effective Hamiltonian 
function TTeff of the dynamical variables a and R. 

Once we have the Hamiltonian, we may calculate the statistical average of a dy¬ 
namical variable of the type F{a, R) simply by 

J dpRdpadadRF{a, 

/ dpn dpa da dR 




However, because the effective Lagrangian is of the form 





R)K 



(26) 


where K = K[a, R) is a symmetric matrix, after the momentum integration the aver¬ 
age can be cast into the form 

J dadRF(a,R)Vd^e-f^^ 

a,R)) = - , - 

JdadRVd^e-^^ 




The matrix 
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and the pseudopotential V are given separately for the two regions. For 0 G [a_, a+] 
we obtain 
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and 


Rc{o) _y _ 2a^A^RA\ + 2a^AB?mJ'i{mR) + -a^R^wAl{2mR) 


dvr 


+ -m^a^A^R^A^ + rAa^AR'^B^ H— rn^a^R^ - aT 

2 2 2 6 3 Rc{a) 


A^R^Al 


—aT 


Rc{a) 


1 A 
+-A- 


4 i?c(a)2 


A'^R^B^ - aT 


A^R^Ai + A 


Rc{a) 


AR^Bl - -aT-^R^ 


9 i?c(a) 


A^R^B^ + -A 


i?c(a)' 


A'^R^B^ 


2 i?e(a)2 


+A- 




(30) 


i?c(a)2 12 Rc{a)^' 

Note that in Eq. (|30D the a -dependence of Rc has to be used explicitly because the 
critical behaviour is determined from it. For the region where (j) ^ [a_, a+] the 
corresponding functions are given by 

An = 


Ki 2 = aR?’mI{mR) 
K 22 = R^ rA I {2mR) 


(31) 


and 


—V = -a^R^m^I{2mR) + -tAo^R^ - -aTa^R^ + —XAR^. (32) 

dvr 2 ^ ^ 6 9 12 ^ ’ 


A number of shorthand notations have been introduced in the previous equations: 

A = R^{a) - R (33) 
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The range of integration for R posses un upper limit given by thermalization. Mo¬ 
tivated by the fact that thermal fluctuations can generate conhgurations with spatial 
size comparable to the critical bubble radius, which may affect the dynamics of a hrst 
order phase transition, the authors of ref. |]T] have estimated the lifetime of fluctua¬ 
tions of an on-shell Higgs held with zero momentum (po = P = 0). This choice 

rehects the fact that critical bubbles are typically much larger than the interparticle 
distance ~ 1/T in plasma. Writing po = oo — i 7/2, one hnds that the dispersion 
relation is 

,2 , , 1.2 
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PI 


+ m2(T) + -y, 


where 
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T^^) being the two-point function for the Higgs held. 

The imaginary part arises at one loop level, but because of kinematical constraints, 
the two loop contribution is actually dominant in the region of physical couplings. The 
thermalization rate 7 for small amplitude scalar huctuations and large spatial size, 
R ~ |p|~^ 7 ~^, is estimated to be of the order 7 ~ 10“^ T near the critical 

temperature, i.e. much larger than the typical hrst order transition time. This means 
that all small amplitude huctuations with size larger than 

7?max = 0 ( 1 / 7 ) (40) 


will ehectively be absent from the mixture of subcritical bubbles and must not counted 
in the thermal averages. In practise, the limit Eq. (^OD is of the order of few times Rc, 
depending on the actual value of 7 . Even if it is not precisely known, its inclusion in 
the calculations is important. Without it all statistical averages would be dominated 
by inhnite, inhnitesimally small huctuations. Technically this can be seen from the 
Eq. (p7|), where the integrals diverge in the limit a —> 0, i? —> 00 . It is important to 
note that the divergence is not a problem of our Ansatz but merely a more general 
phenomenon, which seems to be related to the general infra-red instability problems 
emerging in the calculations of the ehective action. 

We have computed the average initial radius and the amplitude of huctuations at 
T = Tf from Eq. (pT]) numerically using a cut-oh Rmax — 3.3 Rc (we have checked 
numerically that results do not change signihcantly for diherent choices of i?max)- 

We have taken a = 0.048 and varied A between 4 x 10“^ and 10“^. For larger 
values of A the hrst order electroweak phase transition is close to a second order and 
does not proceed by critical bubble formation. 







For instance, when onr phenomenological potential Eq. m is fitted to the two 
loop resnlt for the effective potential calcnlated in [^] for the Higgs mass Mh = 70 
GeV, this yields A ~ 0.061. One can readily verify that the thin wall approximation 
is valid in the chosen range for A. 

In Fig. 1 we show the ratio between and the inflection point 0inf(Tj) as a 

fnnction of A. This ratio is always of order of 0.5. We have also compnted nnmerically 
the ratio {a)/0inf(Tj) which tnrns ont to he in the range (0.13-0.17) in the given range 
for A. We recall that the oscillon stage can be present only if snbcritical bnbbles are 
thermally nncleated with initial amplitnde above the inflection point @]. 

In Fig. 2 we show in nnits of Rc and we have nnmerically compnted (i?) to 

be in the range (1.74-1.58) Rc- Oscillons can be formed only if > i?min — 3/'^(7/) — 
(0.55 - 0.27) Rc i. 

Note that, when A increases, the average amplitnde at the core also increses, 
whereas the smallest available amplitnde to settle into an oscillon stage decreaes 
[§]. However, in spite of this tendency, the avarage amplitnde is always smaller than 

for any chosen valne of A. 

Thns, even if onr resnlts seem to indicate that condition ii) is satisfied, i.e. snbcrit¬ 
ical bnbbles can be thermally nncleated with snfficiently large average initial radins R 
to give rise to the oscillon regime, nevertheless initial average amplitndes at the core 
do not satisfy condition i) since they always resnlt to be smaller than the inflection 
point 0inf(E/). 

Thermal flnctnations are certainly present at the onset of the electroweak phase 
transition, bnt the most probable snbcritical confignration generated aronnd the criti¬ 
cal temperatnre, even if with snfficiently large size, is charaterized by an amplitnde too 
small to begin the oscillon stage. From this resnlt we can infer that the dynamics of 
a weak electroweak first order phase transition is not affected by the presence of long- 
lived oscillons, snggesting that the electroweak baryogenesis scenarios are still viable 
to explain the generation of the baryon asymmetry in the early Universe. However, 
we feel that an important issne deserves fnrther stndy: relaxation time-scales of snb¬ 
critical conhgnrations depend on the natnre of the stochastic force and the strength of 
dissipation provided by the snrronnding thermal bath and their complete knowledge is 
needed to decide if degrees of freedom are in eqnilibrinm or not inside the snbcritical 
bnbble. 

One may realize that this is a crncial qnestion by reminding that the effective 
potential ([l0|) , used to describe the free energy associated to the fluctuations, is usually 
obtained integrating out fermionic and the bosonic degrees of freedom of the theory. 

In performing such a calculation, it is commonly assumed that fermions and bosons 
do have equilibrium distributions with a 0(x, t) background dependent mass. This is 
true only if their interaction times with the background 0(x, t) are much smaller than 




the typical lifetime of the subcritical bubble. 

Since this condition is not always satisfied, a full non-equilibrium approach is 
needed. The latter, however, seems to conhrm, or even strengthen, the results of 


this paper about oscillons [12 
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Figure captions 

Figure 1 The plot of the ratio (a^)^/^/0inf(T/) as a function of A and for a = 0.048. 
Figure 2 The plot of in units uf i?c as a function of A and for a = 0.048. 



